Photonic quantum transport in a nonlinear optical fiber 
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We theoretically study the transmission of few-photon quantum fields through a strongly nonlinear 
optical medium. We develop a general approach to investigate non-equilibrium quantum transport 
of bosonic fields through a finite-size nonlinear medium and apply it to a recently demonstrated 
experimental system where cold atoms are loaded in a hollow- core optical fiber. We show that 
when the interaction between photons is effectively repulsive, the system acts as a single-photon 
switch. In the case of attractive interaction, the system can exhibit either anti-bunching or bunching, 
associated with the resonant excitation of bound states of photons by the input field. These effects 
can be observed by probing statistics of photons transmitted through the nonlinear fiber. 
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Quantum dynamics of strongly correlated systems far 
from equilibrium is a new frontier of many-body physics. 
Intriguing phenomena involving such dynamics have re- 
cently been observed in diverse physical systems ranging 
from ultra-cold atoms [H, 13, 0, S] to individual spins in 
semiconductors @, At the same time, recent experi- 
ments using ultra-cold atoms and optical photons 0, Q 
have opened the door for studies of a novel form of quan- 
tum transport involving strongly correlated photons. It 
has been predicted that these systems can allow for sur- 
prising behavior, such as the dynamical creation of a 
"crystal" of photons [9]. While many similarities exist 
between these photonic systems and condensed matter 
systems involving massive particles, the photonic systems 
also present a unique set of challenges. In particular, 
they do not thermalize, and are inherently open, driven 
systems, which highlights the need to develop novel tech- 
niques for analysis. 

In this Letter, we describe a general technique to study 
the quantum transport of a few field quanta through a 
finite-length, strongly nonlinear one-dimensional waveg- 
uide. This technique allows us to determine the full 
spatial wavefunctions of the photons inside the waveg- 
uide as well as correlation functions of the outgoing re- 
flected and transmitted light. We consider an optical 
waveguide in which the tight confinement of photons near 
the diffraction limit 0, [ll| and the large number 
of atoms with which they interact should enable large 
optical nonlinearities at the single-photon level, which 
is necessary for applications like single-photon switching 
and photonic quantum gates [l2|, [isl, [l^. Unlike non- 
linear optical effects in cavity quantum electrodynam- 
ics (QED) [15], where only a single spatial cavity mode 
is involved, these waveguide systems are more difficult 
to treat in that they contain a large number of spatial 
degrees of freedom, much like low-dimensional, strongly 
interacting condensed matter systems [16|, 
containing a few quanta. As a specific application, we use 



FIG. 1: Attractive {k < 0) and repulsive {k, > 0) photons: 
Output correlation function g2{r = 0) as a function of non- 
linearity. The two-photon wavefunctions {\(l){zi, Z2)\) for four 
values of the nonlinearity are shown. The system size is d=30. 



the calculated reflection and transmission amplitudes to 
demonstrate how such a system can be used to realize a 
single-photon switch. Our analysis reveals the tendency 
for photons to "organize" themselves due to an effective 
tunable repulsive or attractive interaction [as shown in 
Fig. [1]. In the latter case, we show that two-photon 
bound states can form inside the waveguide, which can be 
observed using standard quantum optical measurements 
of the output light. 

Specifically, we are interested in the situation where 
the dynamics of photons inside the waveguide is governed 
by the quantum nonlinear Schrodinger equation (NLSE), 



1 d^^ 

2m dz'^ 



(1) 



Here ^ is a photonic field annihilation operator whose 
specific form will be specified later. The presented frame- 
work can be used to treat the quantum transport of any 
NLSE system (such as interacting atomsjthat are con- 
strained to move along one dimension [l8|, Ell), a 
problem that in general remains relatively unexplored. 
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However, we focus on a photonic system implemented re- 
cently in Ref. where an ensemble of atoms is loaded 
in a hollow-core fiber to mediate the interaction between 
photons. Moreover, we take advantage of quantum op- 
tical techniques [22|, [13] based on electromagnetically in- 
duced transparency (EIT) and photon trapping which al- 
lows us to dynamically control the effective mass and in- 
teraction strength in Eq. ([T]) [9] . The light inside the fiber 
interacts with an ensemble of atoms (with linear density 
no and length L) whose level scheme is shown in Fig. O 
Following Refs. [1, H^, we use forward- and backward- 
going dark-state polariton operators which describe 
excitations of slowly- varying photonic fields £± that are 
coupled to an atomic spin wave excitation cjac through 
a standing- wave control field 0^±{t) = The stand- 
ing wave creates an effective Bragg grating that couples 



^± together and traps them inside the medium [23|, 
much like an optical cavity. In the limit of large opti- 
cal depth (OD = noL^^), the symmetric combination 



(- 



) becomes the only independent quantity 
while the anti-symmetric combination A = ( ^+^^~ ) adi- 

abatically follows, A —^dz"^. The detuning Ai leads 
to a quadratic dispersion that can be interpreted as an 
effective mass m = |(1 + ^ 2\Ai\ )• presence of the 

additional atomic state \d) leads to an optical nonlinear- 
ity [25] whose strength is given by k = 4^xf+?r72)' where 
TiD is the spontaneous emission rate into the guided 
modes and F is the total spontaneous emission rate of 
both states \b) and \d) (Fii)<F). We are primarily in- 
terested in the limit |Ai^2|^F such that m, are mostly 
real and the evolution is dispersive. Finally, we note that 
Eq. ([T]) is written in dimensionless units; all lengths and 
time are scaled by Lcoh = 2(A^ + F^/4)/Fii)no| Ai| and 
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FIG. 2: (a) Four- level atomic system yields a field evolution 
equation in the form of an NLSE. (b) A schematic of quan- 
tum transport inside a finite-length waveguide. An input field 
^+(0) is injected into one end of the waveguide. Inside the 
medium, the effective Bragg grating couples the forward- and 
backward-propagating fields together and the evolution of the 
coupled fields is described by the NLSE. A field is transmit- 
ted (reflected) at the end of the waveguide z = d {z = 0). 
(c) Linear transmission spectrum — 0) as a function of 
detuning, for a system of optical depth OD=noL^^ = 1000. 



tcoh = I Ai 1/211^, respectively. 

First, we investigative the linear transmission {k, = 
0, ^ = (^)) spectrum of a finite system of dimensionless 
length d = 77^, where an input field ^+(0) is driving one 
end of the system [see Fig. [2fb)l . By solving a set of cou- 
pled mode equations for ^± [2J], one can show that like 
an optical cavity, the coupling of ^± inside the waveguide 
creates a set of transmission resonances in the output 
field ^+(<i), which due to the quadratic dispersion of the 
system occurs at 6 = (^)^ [26]. Here 6 = Astcoh is the 
dimensionless two-photon detuning and n is a positive 
integer [Fig. [SJc)]. For a fixed resonant optical density 
(OD d^^), the system losses can be tuned to yield 
either a sharp, low transmission peak (large d) or a high 
broad one (small d) as shown in Fig|2l^c). The finesse of 
this effective cavity is proportional to OD [26] . For a clas- 
sical field, addition of the nonlinear term shifts the trans- 
mission peaks in frequency in an intensity-dependent way 
to the left or right depending on the sign of the nonlinear- 
ity coefficient. Transport in the classical NLSE system 
has already been extensively studied, where the nonlin- 
earity is sufficiently weak that the operators in Eq. [1] 
can be replaced by complex numbers representing their 
mean-field values. In the context of atomic Bose-Einstein 
condensate transport 13, 0, |2^, this yields the well- 



known Gross-Pitaveskii equation, while in nonlinear op- 
tical fibers the classical NLSE explains, for instance, the 
formation of optical solitons [30] . While the nonlinearity 
in our system can be made small enough to also observe 
such effects, here we are primarily interested in explor- 
ing the novel regime where the nonlinearity becomes ex- 
tremely strong. In particular, when a^OD ^ 1, the fre- 
quency shift is predicted to be significant at intensities 
corresponding to a single-photon level. In this regime, 
the quantum transport of a few photons is expected to 
behave fundamentally differently than predicted by clas- 
sical calculations, as we describe below. 

In order to study quantum transport properties of the 
system, we investigate the evolution of the wavefunction 
of the field in the Schrodinger picture. We are interested 
in the situation where the system contains no more than a 
few photons, which is sufficient to describe the properties 
of a single-photon switch. This allows us to truncate 
the Hilbert space so that only vacuum, single-photon, 
and two-photon states are present (the analysis can be 
extended to truncating at any small number of photons). 
We thus write the general state as: 



= I dZidZ2Hzi,Z2,t)^H^l)'^H^2m 

dzO{z,t)^'^{z)\0)^e\0). (2) 



The first, second and third terms on the right correspond 
to the two-photon, single-photon and vacuum compo- 
nents, respectively. Importantly, directly solving for the 
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two-photon wavefunction (l){zi, ^2, t) ahows us to charac- 
terize any non-trivial spatial order between two photons 
in the system. The evolution of the two-photon wave- 
function (j){zi^ Z2^t) under Eq. ([T]) is given by 

idtct> = [dl^ + O + 2k#(^i - Z2), (3) 

while the single-photon state evolves as a free, massive 
particle. Generally, direct integration of such evolution 
equations should be a viable approach for any transport 
problem involving only a few photons (or any other par- 
ticles). A novel characteristic of photonic systems, how- 
ever, is that light inherently can also enter and leave. 
Thus a remaining task is to formulate an "input-output" 
formalism that allows one to relate the wavefunctions in- 
side the waveguide {^<z<d) to the fields injected into 
and propagating out from the system and any of their 
correlation functions, which we describe below. 

Since the NLSE Hamiltonian commutes with the field 
number operator (^"f"^), manifolds with different field 
quanta are decoupled from each other inside the medium. 
However, the system is driven with an input field at 
z = 0, and different manifolds can be coupled at the 
boundaries. In particular, for a classical input field, 
the boundary conditions correspond to that of a co- 
herent state, ^+(z = O)IV^(t)) = aQe-'^^\^l)(t)) where 
ao is the coherent state amplitude and 5 the detun- 
ing. Since ^+ annihilates a photon, the boundary con- 
dition relates any manifold with n photons to that with 
n — 1. Using the fact that the anti-symmetric combi- 
nation A adiabatically follows ^, we can write ^± = 
-1=(^ ^ idz"^) (here we consider the mass to be real). 
Therefore, using the multi-mode wavefunction of the sys- 
tem from Eq. ([2]), the boundary condition for z = 
reads (/>(0, ^2, t) - iS,,(/)(zi, ^2, t)|,,^o = ^e-'^^0{z2,t) 
and 0{z = 0, t) - idzO{z = 0, t) = V2aoe-'^^e, and e 1. 
Similar boundary conditions apply at 2: = d with an in- 
put ao = 0. The intensity or any other field correlation 
function can easily be obtained from the photon wave- 
functions. For example, the normalized second-order cor- 
relation function of the output field ^+(^ = d) is given 

]q ^ (q\ ^ 4|0(d,d)p 

^ ^ (l^(rf)|H4 / dzmz,d)\'Y' 

Before studying the dynamics, we analyze the funda- 
mental modes and energy spectrum supported by the sys- 
tem in the absence of any input field {ao = 0), which can 
be found via the Bethe ansatz technique [3, [3l| . This 
ansatz specifies that the eigenstates consist of a superpo- 
sition of states in which colliding particles exchange their 
wavenumbers k^. Unlike the typical formulation, the val- 
ues of ki here can be complex to refiect the open nature of 
our boundary conditions, which allow particles to freely 
enter or leave. Imposing the boundary conditions, one 
finds that the fundamental modes satisfy: 

^2ik,d ^ (fej + 1)^ {kj - kj + iK.){ki + kj + m) 
(/c^ 1)^ (^i kj ZAc) (/c^ ~\~ kj zhh^ 



where (i, j) can be (1,2) and i j (for detailed deriva- 
tion and generalization for many-body case see Ref. ^). 
Correspondingly, the energy of each state E = /c^+Zcl can 
acquire an imaginary part describing its rate of exit. One 
class of allowed solutions corresponds to the case where 
each photon is freely propagating. In the strong interac- 
tion limit, these wavevectors exactly match the allowed 
values for the non-interacting case, regardless of whether 
the interaction is attractive {hz > 0) or repulsive {tz < 0). 
In the strongly repulsive case, this phenomenon is well- 
known as the fermionization of hard-core bosons [21I] . For 
a large system {d ^ 1), the non-interacting wavevec- 
tors and energy can be approximated by (^, ^) and 

(^)^ + (^)^, respectively. A second class of solutions 
appears for negative k and corresponds to bound states 
whose momenta for large system size can be approxi- 
mated by (/ci, /C2) — ^ ± . In other words, the inter- 
action Hamiltonian V = —2\hz\6{zi — Z2) admits a single 
bound state in the relative coordinate of the two pho- 
tons, while the allowed momenta in the center-of-mass 
coordinates are quantized and determined by the system 
size. This yields a discrete set of bound state energies 

We now examine the properties of the light transmit- 
ted through the waveguide, first considering the repul- 
sive regime n > {). We fix the input field detuning to 
^res = (7r/<i)^, which corresponds to the first transmis- 
sion resonance in the linear regime. Since we wish to 
study an open, driven system out of equilibrium, con- 
ventional analytical methods such as Bethe ansatz or 
quantum inverse scattering [13] cannot be applied here, 
whereas Eq. ([3|) with proper boundary conditions can 
easily be numerically solved. Although all numerical re- 
sults correspond to a specific set of parameters (system 
size, detuning, etc.), the conclusions are quite general. 
In the presence of linear absorption, only a fraction of 
the single-photon state is transmitted and the rest is 




FIG. 3: (a) Repulsive photons: Correlation function g2{r — 0) 
of the transmitted light as a function of optical density (OD = 
noL^^) and single- atom cooperativity when the frequency is 
set to the single-photon transmission resonance with T :^ 20% 
and ^ = 5. (b) Attractive photons: Bound states energies 
(solid) become resonant with energies of incoming photons 
(dashed) for specific values of The error bars represent the 
decay rate of each bound state. 
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dissipated (T < 1); this results in the imperfect trans- 
mission shown in Fig. [2](c), since the single-photon state 
obeys linear dynamics. We note that since the system 
is weakly driven (aod)'^ ^ 1, the noise effect due to ab- 
sorption can be ignored [26|. On the other hand, the 
two-photon state transmission is further suppressed due 
to an extra nonlinear phase shift that brings this state 
out of transmission resonance (for a similar effect in a 
cavity, see Ref. [HI), which results in anti-bunching of 
the transmitted light (^2(0) < 1). The anti-bunching 
implies that the transmitted field acquires non-classical 
character [32[, and the effect becomes more pronounced 
as OD is increased, which increases the effective system fi- 
nesse [Fig.[3l^a)]. Physically, this behavior occurs because 
transmission of the single-photon component of the field 
remains intact while two-photon and higher components 
become strongly reflected, and thus it can be said that 
this system operates as a single-photon switch. More- 
over, we observe a significant spatial deformation of the 
two-photon wave function inside the waveguide as is 
increased, as seen in Fig. [H In particular, a cusp de- 
velops along the diagonal zi = Z2 and the majority of 
the photon density lies in the off-diagonal regions, in- 
dicating that the two photons repel each other. In the 
limit of hz^oo^ the diagonal elements must be completely 
suppressed which leads to perfect anti-bunching. Similar 
behavior involving the "self-organization" ofphotons in 
an NLSE system has been discussed in Ref. [y]. 

In the case of attractive interaction, the system ex- 
hibits a very different behavior from the repulsive case. 
In Fig. [11 we plot ^2(0) for the transmitted field versus 
Kd. For small \hz\d^ the nonlinear phase acquired by multi- 
photon components enhances (suppresses) their reflec- 
tion (transmission), much like in the repulsive case. This 
is responsible for 5^2(0) dropping slightly below unity, de- 
spite the attractive interaction, and yields a weak single- 
photon switching effect. At larger values of \f<i\d^ oscilla- 
tions develop in the correlation function, yielding strong 
bunching behavior at particular values of tvd. A closer 
analysis reveals that the locations of these resonances 
correspond to resonant excitation of specific two-photon 
bound states, 2Sres = E^- In particular, while the de- 
tuning 5res of the input field is fixed, the bound state 
energies vary quadratically with changing a^, causing 
the two to come into resonance at particular values of 
K. This is shown in Fig. [3](b), where the bound-state 
energies are evaluated via Eq. This effect is further 
confirmed by examining the two-photon wavefunction at 
each of these oscillation peaks [Fig.[T], where the photon 
density becomes localized along the diagonal and indi- 
cates a bound state in the relative coordinates. Simul- 
taneously, an increasing number of nodes and anti-nodes 
develop along the diagonal for increasing \i<i\d^ which are 
associated with the higher momenta of the center-of-mass 
motion. These resonances deviate significantly from the 
predictions of the semiclassical picture and also have no 



counterpart in cavity QED, due to the unique spatial de- 
grees of freedom available in our system. Importantly for 
experiments, these bound states can be probed with clas- 
sical light, simply by examining higher-order correlation 
functions in the output field. For example, the bunching 
peak associated with the second bound state \n\d c^b^ 
can be observed with parameter choices T 1% and 
^ = —5, using an optical density of OD 3500 and 
coupling efficiency =0.2. 

In conclusion, we have presented a theoretical ap- 
proach to investigate transport properties of a few field 
quanta inside a one-dimensional, finite-size nonlinear 
medium obeying the NLSE. When the system is driven 
out of equilibrium with a coherent field at the boundary, 
the strong interaction inside the medium is manifested in 
the correlation functions of the transmitted field. In par- 
ticular, for a repulsive interaction, the transmitted field 
is anti-bunched, while for an effective attractive interac- 
tion, the transmitted field can be either bunched or anti- 
bunched. Moreover, in the attractive case, we show that 
the bunching is due to the resonant excitation of photonic 
bound states. Recent progress in loading a large number 
of cold atoms into hollow-core optical fibers [7] should 
allow one to probe these effects in ongoing experiments. 
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